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The Bipolar Spherical Harmonics (BipoSH) form a natural basis to study the CMB two point 
correlation function in a non-statistically isotropic (non-SI) universe. The coefficients of expansion 
in this basis are a generalization of the well known CMB angular power spectrum and contain 
complete information of the statistical properties of a non-SI but Gaussian random CMB sky. We 
use these coefficients to describe the weak lensing of CMB photons in a non-SI universe. Finally we 
show that the results reduce to the standard weak lensing results in the isotropic limit. 


I. INTRODUCTION 

We assume the CMB temperature anisotropies to be 
Gaussian, which is in good agreement with current CMB 
observations. Consequently the two point correlation 
function contains complete information about the un¬ 
derlying CMB temperature field. Generically, the two 
point correlation function can be expressed in terms of 
the spherical harmonic coefficients of CMB temperature 
maps as follows, 

C(ni, n 2 ) = (AT(ni)AT(n 2 )), 

= X (aima^YUn^Y^,^), (1) 

Iml'm' 

where (a 2 m a*, m ,) is commonly referred to as the covari¬ 
ance matrix. 

Now, if one assumes the universe to be statistically 
isotropic (SI), it can be argued that the two point cor¬ 
relation function of the CMB temperature field cannot 
have any explicit dependence on the directions fii and 
n 2 and hence can only depend on the angular separation 
between the two directions, 

c(h 1: h 2 ) = c(h 1 -h 2 ) = c(0), ( 2 ) 

where 0 = cos~ 1 (fii ■ n 2 ) is the angular distance between 
the two directions. This correlation function can be ex¬ 
panded in the Legendre polynomial P; basis as follows, 

9/ i i 

C^fii, n 2 ) = X 47r Cf-FHni • n 2 ), (3) 

i 

where CMB angular power spectrum Ci is the coefficient 
of expansion in this basis. It can be shown that Ci is 
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related to the covariance matrix through the following 
expression, 

Cl — ((LlmQ'l 1 m' )^ 22 ' &mm' i ( 4 ) 

implying that in a SI universe the complete statistical 
information on the CMB field is encoded in the diagonal 
elements of the covariance matrix. 

II. DESCRIBING A NON-STATISTICALLY 
ISOTROPIC UNIVERSE 

If we make no assumption about a SI universe, then 
the two point correlation function explicitly depends on 
the directions ni and n 2 . This property reflects itself in 
harmonic space, as non-vanishing off-diagonal elements 
in the covariance matrix. This implies that the angu¬ 
lar power spectrum by itself cannot fully characterizes 
the statistical properties of a non-SI CMB sky. In fact, 
searches for non-vanishing off-diagonal elements in the 
covariance matrix forms the basis of all tests searching 
for violations of isotropy. 

In a non-SI universe the two point correlation is best 
expressed in terms of a basis which has this explicit bi¬ 
directional dependence. The natural choice is then that 
of the Bipolar Spherical Harmonic (BipoSH) functions [1, 

2] , which form a complete orthonormal basis for functions 
defined on S 2 x S 2 . The two point correlation function 
for any non-SI field can be expanded in the BipoSH basis 
as follows, 

C(ni,n 2 )= X A hhi Y h( A t) ®Yi 2 (n 2 )} LM . (5) 

LMhh 

where Af^ are the coefficients of expansion in the Bi¬ 
poSH basis. The BipoSH basis functions themselves, can 
be expressed in terms of the standard spherical harmonic 
functions Y) m as, 

{Y^n^Y^^LM = X C “il2m 2 ^im 1 (fil)^! 2 m 2 (n 2 ), 

7711771-2 


( 6 ) 
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where are the Clebsch-Gordon coefficients. We 

refer to L as BipoSH multipoles and Zi, Z 2 as spherical 
harmonic (SH) multipoles. The Clebsch-Gordon coef¬ 
ficients are defined only when their indices satisfy the 
following constraint equations, 


|Zi — Z 2 I < L < h + I 2 , 

(7a) 

mi + 7772 = M , 

(7b) 

VI 

VI 

T 

(7c) 

— 12 < 7772 < h ■ 

(7d) 

Using the orthonormality property of the BipoSH 
functions, it can be shown that the BipoSH spectra 
are related to the the covariance matrix as follows, 

basis 

aLM 

‘ 1^2 

\LM _ „ \ r LM 

^hh ~ 2_^ \ a hm 1 ai2m 2 /L'l 1 m 1 l2m2 ' 

(8) 

m\m 2 


As evident from the above equation, all the information 
in the covariance matrix is captured in the BipoSH spec¬ 
tra. Specifically it can be shown that the BipoSH spec¬ 
trum A™ is related to the CMB angular power spectrum 

as follows, 


A™ = (-1 )%C l , 

(9) 


III. WEAK LENSING OF THE CMB 
TEMPERATURE FIELD 

The CMB photons reach the observer after travers¬ 
ing through the intervening large-scale structures (LSS). 
This results in the CMB photons deviating from 
their geodesics in the unperturbed Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric. We work with the 
approximation that the deviations from the geodesics are 
small, such that the photons sample the same gravita¬ 
tional potentials as they would have along the unper¬ 
turbed trajectories, termed the Born approximation [ 8 ]. 

Weak lensing of the CMB photons results in the CMB 
temperature map getting remapped, 

T(n) = T(n + A), (12) 

where A is the deflection to the direction of photon ar¬ 
rival. On solving the photon geodesic equation under the 
Born approximation it can be show that A is related to 
the projected lensing potential ip as follows, 

A(n) = VV’(n) • (13) 

The projected lensing potential ip is expressed as a 
weighted sum of the gravitational potential encountered 
by the photon along its trajectory, 


where n; = V2 1 + 1. Therefore, the BipoSH spectra are 
a generalization to the commonly studied CMB angu¬ 
lar power spectrum. To summarize, while the BipoSH 
spectrum A^ characterize the isotropic component of the 
field, the rest of the BipoSH spectra Ajpj^( L / 0), de¬ 
scribe the non-SI component of the field. The primary 
scheme to detect a violation of isotropy in the data has 
been to search for these non-vanishing BipoSH spectra 
with L > 0 [3-6]. 

Before we proceed, we discuss a few more interesting 
properties of the BipoSH spectra. For any fixed Zi, Z 2 
the BipoSH spectra A LAI have the same mathematical 
properties as the spherical harmonic coefficients of a map. 
Motivated by this property it is possible to define the 
reduced BipoSH map formed by summing over the SH 
multipoles of the BipoSH spectra [7], 

n LM = Y,Kh- ( 10 ) 

hh 


It is also possible to define the Bipolar power spectrum 
(BiPS), a coordinate independent quantity which is a bi¬ 
ased estimate of SI [1] and is defined as 

«L= E \ A hh I' ■ (11) 

h ,h,M 

We will discuss how these completely generically defined 
quantities are related to some cosmological observables 
in a later section of this article. 


V’(n) 



vo-v 

. VoV . 


0(»7,n), 


(14) 


where 77 is the angular diameter distance which is identi¬ 
cal to the comoving distance in a flat universe and <^( 77 , n) 
is the gravitational potential encountered by the photon 
coming in from direction n at a comoving distance of 77 . 

Analogous to the CMB temperature field the projected 
lensing potential is a Gaussian random field defined on a 
sphere and thus in harmonic space 


ip(h) 

= EE fcnYfcn ( fi ) ’ 

(15) 


kn 



— ('iftkn'lftk'n'^^kk' $nn' i 

(16) 


where ipkn denotes the spherical harmonic coefficients of 
the lensing potential and C^ the corresponding angular 
power spectrum in a SI universe. 

In a non-SI universe the statistical properties of ip are 
characterized by the following BipoSH spectra, 

Kl = E (V’UmV’W « lfe2 „ 2 . (17) 

n 1 n 2 

where the BipoSH coefficient (V K ^ 0) encodes 

the information in the off-diagonal components of the 
harmonic space covariance matrix of the projected lens¬ 
ing field. 

Weak lensing of the CMB photons through the inter¬ 
vening LSS leave measurable imprints on the CMB two 
point correlation function. Specifically lensing introduces 
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coupling between different multipole moments of the tem¬ 
perature field, which are absent in an unlensed SI CMB 
sky (See Eq. 4). The lensing modification to the CMB 
angular power spectra (the diagonal components of the 
covariance matrix) have been well studied [9, 10]. 

In the following sections of this article we discuss weak 
lensing induced modifications to the BipoSH spectra. 

IV. BIPOSH REPRESENTATION OF A LENSED 
CMB SKY 

Since we measure the lensed CMB sky, while modeling 
the two point correlation we must account for the weak 
lensing induced modifications to the statistics of the map. 
As discussed in the previous section, lensing remaps the 
temperature field. Given the statistics of the projected 
lensing potential ip, it is possible to evaluate the lensing 
modifications to the BipoSH spectra. 

Since the deflection A is small, we can Taylor expand 


and express the lensed temperature field T(n) in terms 
of the unlensed temperature field T(n) and the deflection 
field A (n). In order to evaluate the lensing modifications 
to the two point correlation of the CMB temperature field 
we retain terms which are at most second order in deflec¬ 
tions. Hence to leading order terms in the deflection field 
A (rc), the lensed temperature field T can be expressed in 
terms of the unlensed temperature held T as follows, 

T(n) = T(n + A), 

« T(n) + V°V>(fi)V 0 T(n) 

+ iv a ^(n)vV(fi)V a V 6 T(n) + 0(V> 3 ). (18) 

In the above equation we have used the fact that the 
deflection held A(n) is related to the gradient of the pro¬ 
jected lensing potential. Consequently, the two point cor¬ 
relation of the lensed temperature held can be expressed 
as follows, 


(T(n!)f(n 2 )) = (T(n 1 + A 1 )T(h 2 + A 2 )> 

= <T(n 1 )T(n 2 )) + (V“^(fii)V„T(n 1 )T(n 2 )) + (V“^(n 2 )V a T(n 2 )T(n 1 )) 

+ (V a V’(fii)V 6 V'(fi 2 )V Q T(n 1 )V b T(n 2 )) 

+ i(V a V’(fii)VV(fii)V a V 6 T(n 1 )T(n 2 )) + i(V a V’(n 2 )V^(n 2 )V a V 6 T(n 2 )T(n 1 )) . (19) 


Here we follow the Einstein summation convention, i.e., 
repeated indices of sky coordinates are summed over. 

We make a small diversion to discuss some subtle as¬ 
sumptions we have made in the analysis that follows. 

• We emphasize the point that unlike the CMB tem¬ 
perature held, we do not treat the projected lensing 
potential as a stochastic held, instead we treat it as 
a constant held around an observer. This is justi¬ 
fied by the fact that we are surrounded by a par¬ 
ticular realization of the LSS and we see the CMB 
through this particular realization. 

• Since the averaging is done over a single realization, 
the term linear in the projected lensing potential ip 
does not vanish. 


The remarks made above are very critical to the anal¬ 
ysis presented in this article. Since the dominant con¬ 
tribution to the projected lensing potential comes from 
very local structures, any anomaly in the local distribu¬ 
tion of structures around us would immediately rehect in 
the CMB observations. We demonstrate in the following 
sections that the BipoSH coefficients are observables con¬ 
structed from observed CMB maps which can in principle 
allow us to test the nature of the local universe. 

The operation of taking the gradient transverse to the 
line-of-sight (V a = V?) and the averaging (...) operation 
commute and we use this property to express the lensed 
two point correlation function of the temperature field as 
follows, 


J 


(r(n 1 )T(n 2 )) = 


+ 


(T(n!)T(n 2 )) + V“V>(fii)V" 1 (T(n 1 )T(fi 2 )) + V^^V^T^T^)) 
V| 1 V^(^(n 1 )^(n 2 ))V" 1 Vf(T(n 1 )T(n 2 )) + 1 [V| x (^(n 1 )^(fii))V^ 1 V" 1 

V5 2 V^(^(fi2)^(fi2))V^V" 2 (T(n 2 )T(n 1 ))' . 


(T(n 1 )T(n 2 )) 

( 20 ) 


where the gradient now operates on the ensemble aver¬ 


ages of the field, which can be expressed in the BipoSH 
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basis. Note that the terms which are linear in ip are 
written explicitly without the ensemble averaging. We 
reiterate that while dealing with terms which are first 
order in ip the projected lensing potential is not treated 
as a stochastic field, instead it is thought of as a constant 
field. 

For notational convenience, we identity each of the 
terms in the two point correlations in the above equa¬ 
tion as, 

A = A + (a 1 + a 2 ) + l 3+ (71 + 72 ), 

= A + a + /3 + 7 , ( 21 ) 


where a denotes the terms linear in the lensing potential 
ip, ft denotes terms which are quadratic in ip but con¬ 
sist of only first derivative of the temperature field and 
7 denotes terms which are quadratic in ip but consist of 
double derivative of the temperature field. Next, we ex¬ 
press Eq. 20 in harmonic space, by expressing each two 
point correlation in the BipoSH basis. The transverse 
gradients now operate on these basis functions. A de¬ 
tailed demonstration of this procedure can be found in 
Appendix A. 

Finally this relates the BipoSH coefficients Ap^ of the 
lensed temperature field to the unlensed temperature Bi¬ 
poSH coefficients Ap^ through the BipoSH coefficients 
for the lensing deflection field 


A: 


LM 


= + [ icffi + 2 affi] 


+ [ 


LM 


.XMl 

2 Of p 2 J 


aLM 

7££(* 


LM 
a hl 2 


W’fc' 


m A-i i 


L'M' 


KN A 


L'M' 


LM 


Pi 


■«(* 


k\k 2 ‘ 


A l ' m '\ 

A hU > 

( 22 ) 


where CXpfp is the lensing correction term which is first 
order in ip and (Ppip, 7 are the lensing correction 
terms which are second order in ip. 

In the following sections we evaluate each of these 
terms explicitly and discuss their consequences. 


A. First order correction terms 


We begin by evaluating the lensing correction to the 
BipoSH coefficient arising due to the term which is linear 
in the lensing field ip{n). The linear order correction is 
given by the following expression, 


„,iM , ChF(li,L, l 2 ) + Ci 2 F(l 2 , L, h) L0 

= iPlm -—- v ;— c, 




L'jiO 


n ~hoi 2 o 


, Y^,/, \ ' r LM \ " AL’M'nL'M ' 

‘ / v 'rkn / v L/ /imiZ 2 m 2 / v 

mi m 2 L' M '131713 


kn 


jm\nm3 , aL'M'/iL'M' jm2nm3 
X 1 hkl 3 ' -™-lih ^hm 1 l3m3 1 l 2 kl3 


(23) 


where the first term is assuming the underlying unlensed 
temperature field to be SI i.e. L' ^ 0 while the sec¬ 
ond term is the non-SI term. The integral in 


the above equation can be solved to yield the following 
expression, 

f m i ™ 3 _ / Af., Y. :a Y V* 

1 l 3 kl 3 / UI1 1 1 kn 1 h m 3 :aI hm.! 

_ A(^3, k, l\ ) n fc n ;3 M 1 m 3 (OAi 

\/47r n ; '^kOl 3 Cr'knl 3 m 3 > 

where F{Ia, k, l±) = [^ 3(^3 + 1 ) + k(k + 1 ) — ^ 1(^1 + l)]/ 2 . 

The first thing to note in Eq. 23 is that this term van¬ 
ishes on taking an average over the lensing field ipkn- 
However while making observations of the CMB sky we 
are surrounded by a single realization of the large scale 
structure and hence in practice it is this single realization 
of the LSS that biases our observations. 

In Eq. 23, the first term results from the isotropic com¬ 
ponent of the CMB sky whose statistical properties are 
fully specified by the angular power spectrum. Examin¬ 
ing the first term also yields the fact that the BipoSH 
coefficients derived from the lensed CMB sky, even in a 
statistically isotropic universe, will be rendered non-zero 
due to weak lensing. While the amplitude of the BipoSH 
spectra is determined by the magnitude of the projected 
lensing potential, the spectral shape of the signal is com¬ 
pletely determined by the CMB angular power spectra. 
It can also be seen that there are no corrections to the 
power spectrum arising from this term, since a™ identi¬ 
cally vanishes. 

The second term describes the coupling between vari¬ 
ous BipoSH spectra due to weak lensing. It is not possi¬ 
ble to extract much information from this term without 
assuming a particular model of isotropy violation. 


Estimators for lensing variables. 


As is evident from the above discussion, the BipoSH 
coefficients encode information about the projected lens¬ 
ing potential. It is possible to derive an estimator for 
the projected lensing field and the corresponding power 
spectra in terms of quantities derived from BipoSH coeffi¬ 
cients which can be measured from observations of CMB 
temperature maps. The projected lensing potential is ex¬ 
pressed in terms of the bipolar maps defined in Eq. 10 
through the following expression, 


IpLM = 


\Z4ttIIlTZ 


LM 


[< O h F(h,L , l 2 ) + C h F(l 2 , L, h)] CfU 2 
hh 


<?5) 


where C£, a = H h TL l2 C[f 0l20 . 

Similarly it can be shown that the power spectrum of 
the projected lensing potential is related to the BiPS kl 
defined in Eq. 11, through the following expression, 


_ 


_471-fi^_ 

]T [C h F(h,L, l 2 ) + C l2 F(l 2 , L, h)Cp ll2 } 2 

llh 


.( 26 ) 
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Note that kl is a biased estimate of the lensing power 
spectrum. Hence studying the two point correlation of 
the CMB temperature field in the BipoSH basis very nat¬ 
urally leads to estimators for the lensing variables (the 
lensing field i/jlm and the lensing field power spectrum 
CJ7). These constructions are similar to the lensing esti¬ 
mators proposed in [ 11 ], though the ones presented here 
are not optimal in the sense that they do not minimize 
the variance on the reconstructed lens harmonics i/jlm- 
The above discussion is fully valid for lensing recon¬ 
struction in a SI universe. This scheme of lensing re¬ 
construction is approximate for a non-SI universe. For 
this scheme to be valid in a non-SI universe it is nec¬ 
essary that the second term in Eq. 23 is significantly 
smaller than the first. This may be the case since Ci is 
always larger than the other BipoSH spectra in a non- 
SI universe. While in a SI universe the BipoSH spectra 
measured from the lensed CMB sky can be directly used 
to infer the lensing potential, in a non-SI universe it will 
be necessary to subtract the BipoSH spectra of the un- 


lensed CMB sky from the measured BipoSH spectra in 
order to estimate the lensing potential. 


B. Second order correction terms 

We now evaluate the other correction terms to the 
lensed temperature BipoSH coefficients in Eq. 22 which 
are and Both these terms are second order 

in the projected lensing potential if). In the most general 
scenario, these terms are functions of the BipoSH coeffi¬ 
cients of the projected lensing potential and that of the 
unlensed temperature field. Interestingly, the 7 term, like 
the first order correction term a, comprises of two sym¬ 
metrical terms arising from ensemble averaging along all 
line of sights. However, because of the complexity of the 
second order contribution of the lensing potential, the 
integral in the 7 term is more involved. 


oLM \ ' \XjKN \ ' rtKN \ ' aL'M' \ ' siL'M' \ ' nLM jm\n\m 3 j-m 2 

Pilh / y ^ kxk 2 / y ^'k 1 n 1 k 2 n 2 / , A l 3 li / , U l 3 m 3 Um4 / , '-'hm 1 l 2 m 2 1 l 1 k 1 l 3 1 l 2 h 

KNkik 2 n 1 n 2 L'M'! 3 i 4 m 3 m 4 m 1 m 2 

LM _ LM , LM 

7 hi 2 — 17 hh +2 7qz 2 ! 


7l 2 7Tl4 

k 2 U ’ 


(27a) 


_ v \d kn r KN v c LM 

o / v ^ k\k 2 / j ^ kinik 2 n 2 / v Wimi 


KNk\k 2 


nKN 

■'k\n\k 2 n 2 
n\n 2 m i m 2 


l 2 m 2 


v \ A a L'M' nL'M' jm\n\mzn 2 , a L'M' r>L'M' jm 2 r 

X / v [^3^2 '^hmsl 2 m 2 J l 1 k 1 l 3 k 2 ' ^hh J l 2 ki 


n\rrizn 2 

hk 2 


L' M' 131713 

where the term is the following integral, 




(27b) 


(28) 


The evaluation of this integral is illustrated in Appendix B. This integral can be explicitly evaluated in terms of 
Clebsch-Gordon coefficients and is given by the following expression, 


t? 71 i n 1171,3 n 2 
l\k\l3k 2 


Q'kiQ'k 2 Q'i3 


47T 


n / 


kn 


kn sili—mi 

k\n\k 2 n 2 ^k—nl3—m3 


a h XI fe2 -1 C'olo {! + (-l) fcl+fc2 - fc } + b h X lfe2 -r C l k-21 3 2 {1 + (- 1 ) 


k\-\~k 2 -\-l\-\-l3 1 


(29) 


where a*, = \/k{k + l )/2 and bk = \J(k — 1 )(fc + 2)/2 . 


Reproducing standard lensing results. 


Explicit evaluation of these terms require direction de¬ 
pendent statistics of the primordial fluctuations and that 
of the projected lensing potential which in turn would 
require specific modeling of non-SI. The formalism pre¬ 
sented here is completely generic and hence we refrain 
from discussing specific models of non-SI. 


If one assumes that the unlensed temperature field and 
the projected lensing potential to be SI, then it is ex¬ 
pected that we recover the standard weak lensing mod¬ 
ifications to the CMB angular power spectra. In the SI 
case the only BipoSH modes that would contribute are 
with indices L l = 0 and K = 0 in Equations 27a and 27b 
and we show that this enforces the second order correc- 
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tion terms to reduce to the standard tensing results with 
L = 0 and l\ = / 2 as seen in the following expressions, 


oLM _ 

Pi ih ~ 


X 


^ LM 
Plih 


(- l) h U h S L0 SM 0 S h i 2 J2 c it c h 

kih 


F(h, ki,i 3 ) n fcl n fc3 ^ ii0 ] 
n Zl 1 

(—1) ,1+1 TI^Slo^moSi^ h{h + l)Cii 


(30) 


^1(^1 + l)( 2 &i + F) ^ 
^ 87T kl 


The delta functions in the above equations reduces Eq. 22 
to the SI lensed angular power spectrum for which the 
only non-vanishing BipoSH coefficients are of the form 
j4°i°i = (—1 ) il IIz 1 C'z 1 . The weak lensing correction to 
these terms is only due to the terms which are second or¬ 
der in the lensing potential, i.e. ^ 8 °° and 7 °^. The first 
order correction term does not contribute to the lensed 
angular power spectrum when the unlensed temperature 
field is considered to be SI. Hence using the correction 
terms given in Eq. 30 & Eq. 31 in Eq. 22 yields an ex¬ 
pression for the lensed CMB angular power spectrum Ci 
in terms of the unlensed CMB angular power spectrum 
Ci and the projected lensing potential power spectrum 



C h = C h - Rh(h + 1 )C h + J2 c tt c i 3 

kiI 3 


F(l 3 , ki,h) 
47r 


thJhzrho 


1 2 


n. 


/C1OZ3O 


(32) 


where, 


R = £ h{kl + 1)(2fc i + l) Ct? . (33) 

ki 

This result matches exactly with the expression derived 
for the lensed CMB temperature power spectrum in ex¬ 
isting literature [ 10 ] and serves as a consistency check for 
the formalism discussed in this article. 


V. MODIFICATIONS TO THE CMB ANGULAR 
POWER SPECTRUM 

The most general expression for the lensing modifica¬ 
tion to BipoSH spectra is given in Eq. 22 where the indi¬ 
vidual correction terms are given in Eq. 23, 27a and 27b. 
These lensed BipoSH spectra ■Af'™ arise due to coupling 

of BipoSH spectra A of the non-SI unlensed temper¬ 
ature field and the BipoSH spectra describing the non-SI 
lensing potential. Here we explore the contribution of 
these non-SI terms to the lensed CMB angular power 
spectrum. 

It is generally believed that searching for off-diagonal 
power in the harmonic space covariance matrix is the only 


way to detect deviations from SI. However weak lensing 
mediates a coupling of power between the diagonal and 
off-diagonal components of the covariance matrix. Hence 
in a non-SI universe there could be modifications to the 
CMB angular power spectrum due to mixing of power 
between the diagonal and off-diagonal elements of the 
covariance matrix. In principle these non-standard devi¬ 
ations could be used to test for deviations from SI. 

To quantify the effects of non-SI on the lensed CMB 
temperature power spectrum we evaluate the lensing 
modification to the BipoSH spectrum A]) 0 . Evaluation 
of Eq. 22 for the BipoSH modes L = 0 results in the 
following expression, 


200 
A hh 

C h 




- 4 00 4- 


= c h 


+ 7$ a (*££, A&f) 

(- 1 )' 


n. 


Kl+ztl+A] 


(34) 


We consider three scenarios, first we consider the case 
where the CMB sky at the surface of last scattering is SI 
but the lensing potential is non-SI, next we consider the 
case where the CMB sky at the surface of last scattering 
is non-SI but the lensing potential is SI and finally we 
evaluate the case when both the CMB at the surface of 
last scattering and the lensing potential are non-SI. 

In the first case, the unlensed CMB sky is SI implying 
that only L' = 0 BipoSH spectra are non-vanishing; how¬ 
ever the lensing potential is non-SI which means that all 
BipoSH spectra K £ [0,1, 2,3 • • ■ ] contribute. However, 
using properties of the Clebsch-Gordon coefficients [12] 
it can be shown that setting L' = 0 results in only K = 0 
BipoSH spectra of the lensing potential contributing in 
all the terms a , j3 and 7 . In particular we find that the 
contribution from the first order term a? 0 , vanishes while 

lll2 

the contribution from the second order terms /3i^ 2 and 
7 reduce to the standard lensing terms as in Eq. 32. 

In the next case we have the lensing potential to be SI 
implying contributions only from K = 0 BipoSH spectra 
and a non-SI underlying unlensed CMB power spectrum 
which mean that in principle any BipoSH mode L' £ 
[0,1,2,3-••] could contribute. In this case too, setting 
K = 0 and using the properties of the Clebsch-Gordon 
coefficients, it can be shown that only L' = 0 mode can 
contribute. This again results in the correction terms 
reducing to the standard lensing results. 

Thus we see that the corrections to the power spectrum 
due to non-SI power vanishes on either the lensing po¬ 
tential being SI or the CMB at surface of last scattering 
being SI. These two cases are significant as it shows that 
late time non-SI power generation in the lensing poten¬ 
tial, for example from direction dependent dark energy, 
does not effect the angular power spectrum. Therefore 
any non-standard lensing induced modification to the an¬ 
gular power spectrum in a non-SI universe can only be 
at second order in the parameter characterizing the de¬ 
viation from isotropy. 
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Finally we consider the case where the lensing potential 
as well as the unlensed CMB sky are non-SI. The terms 
in Eq. 34 which contribute to non-standard correction 
to the angular power spectrum Ci are for L' ^ 0 and 


K ^ 0; the other possible values reduce the correction 
terms trivially to the standard lensing results as discussed 
above. Contributions from each of the correction term is 
explicitly evaluated and can be expressed as 


□oo 


v 00 

'lili 


<1 = 2 (-l) ii nq E V’fcA F{h ^l ll) Cofao ■ 

knl 3 V47F 1 fe 

iT ,KN r KN V' aL'M'i i\h+h+lj ^ 3^4 F(h, k 1: l 3 )F(l 1 , k 2 , U) 

^k\k 2 / ; ^k 1 nik 2 n 2 / > ^hU V ^ jt rr~ 

KNk\k 2 nm 2 L'M'hU ll 

nkiO nk 2 0 \ ' /_ -\\m\nh'M' rtl\m\ nl\—m\ 

^'hOhO^'liOhO / j v ^l3m3lAmA^'k\n\lzmzk 2 n 2 l±rn4 L ’ 

mim^rriA 

K^O LVO /_ lVl 

E GyKW \ ^ nKN \ ^ aL'M' \ \ A /_i\mi/iL , M / 

^fci/c 2 / v 0 /einifc 2 n 2 / > jj / > v x / mi 

KNk\k 2 n\n 2 L'M'ls 1 mi 7713 

x J dhY^WY^inW 3^3 :ab (fi)^ imi • 


(35a) 


(35b) 


(35c) 


Thus the non-standard lensing modifications to the 
lensed CMB temperature anisotropy angular power spec¬ 
tra arising due to the universe being non-SI only comes 
about due to the coupling between the off diagonal terms 
of the covariance matrices of the non-SI lensing field and 
the non-SI CMB temperature field (i.e. only through cou¬ 
pling between BipoSH modes with K ^ 0 and L' ^ 0). 


VI. LENSING CORRECTIONS TO COSMIC 
VARIANCE 

In this section, we take a diversion to analyze the ef¬ 
fects of the weak lensing correction terms a, /3 and 7 in 
a SI universe. As argued in section V, weak lensing in a 
non-SI universe can leave imprints of SI violation in the 
CMB angular power spectra (See Eq. 34). However, if 
we were to be sitting in an isotropic universe, then the 
contribution from these non-standard correction terms 
would vanish, i.e. (o^°) = (/3$°) = ( 7 $°) = 0. Note how¬ 
ever that the variance of these terms need not vanish, 
i- e - (K°| 2 ) ^ 0, (I/??/ 3 ! 2 ) ^ 0, <| 7 “I 2 ) 7 ^ 0. These proper¬ 
ties are similar to that of noise, except that this noise is 
of cosmic origin. 

We evaluate the corrections to cosmic variance arising 
because of weak lensing in a SI universe. To recap, the 
weak-lensing contributions to the lensed bipolar coeffi¬ 
cients are 


BipoSH coefficients [13] is given by the following expres¬ 
sion, 


/ a LM a *L'M' \ 
\^l\l 2 1 9 / 


= CiMSwM’, + (-1 ) h+h+L S h v 2 6 l2Vi ] 
x Sll'Smm' + CqC^n^n;/^—l) il+il 
x S^^Si^Slo^moSl'qSm'o • (37) 


The variance of the lensed CMB angular power spectrum 
to first order terms in Cf is given by 


= {C 1 C 1 ) - (Cif (38) 


Since the a™ does not contribute at the power spectrum 
level but has non-zero contribution to the variance, we 
can write the individual terms in Equation 38 as 


(Ci) = Q+^{(AT)+(7ni 

(Cl) 2 = Cf + 2 { — -^(0/3°° + C n °°) + 0[(Cf) 2 ] (39) 

It; 

(Cf) = ({C l + { -^(a™+f3™ + 7D 0 )} 2 ) 

R/ 

= (cf) + 2 ^[(c,<) + (QO + (Caff)} 

lb 

+ iK/°V/°) + o[(^) 2 ]- ( 40 ) 


We evaluate the corrections to cosmic variance arising 
due to the term linear in the projected lensing potential, 
which reduce to the angular power spectrum under SI The variance correction due to this term can be shown to 
Aff = (-lyiliCi. For the case of a Gaussian and SI be linear in the projected lensing potential power spec- 

CMB sky, it can be shown that the covariance of the tra Cf while the correction coming from other terms is 


\LM 

i /,/2 


= A 


LM 

hl2 


LM 

a Ll2 


+0(V> 3 ,V> 4 ), (36) 


..LM 











quadratic in Cf hence is expected to be a sub-dominant. 
As seen in Equation 40, the terms contributing to the 
correction of the cosmic variance and are linear in C'f 
are <|ogf>| a >, ogf>), <A* 00 /?°°> and 

The individual terms are evaluated in Appendix C. Us¬ 
ing these correction terms which are relevant in the case 
of a SI universe, the total cosmic variance is given by the 
following expression, 


Ci 


21 - 1-1 
21 + 1 


- C? 


E°! 


F(l,k,l) 

VAtt 


tt 

LL k^kOlO 


(41) 


A. CMB polarization BipoSH coefficients 

Casting aside the constraint of SI, the most generalized 
CMB polarization two-point correlation function can be 
expanded in the BipoSH basis using spin-2 spherical har¬ 
monics in the following manner, 

C'(n,n') = (±A(n)±X(n / )) 

= <+A(n)+X(n')> + (+A(n)_X(n')) + c.c. 
— ^ ] (E hmi T iB[^m\) 

lm 

x m 2 m 2 + 'i / Bi 2r n 2 ) 2^27712 

H - (Ei 2 m2 ^^ 2 ^ 2 ) —2 ^2^2 ] + c.c., (45) 


VII. EXTENDING THE FORMALISM TO CMB 
POLARIZATION 


where c.c. are the complex conjugate terms. However, 
the two-point correlation function can also be expanded 
in the BipoSH basis as follows 


The higher order weak lensing correction terms dis¬ 
cussed in Section IV will affect the CMB polarization 
field too. As a preliminary study, we lay down the formal¬ 
ism of CMB polarization BipoSH spectra. Weak lensing 
modifications to these coefficients would follow the same 
trend as in the temperature case. 

The CMB polarization field is a spin-2 field and can 
be decomposed into stokes Q(n) and U(n) fields. Us¬ 
ing standard convention of decomposing the polarization 
field into gradient (E) and curl (B) modes which are co¬ 
ordinate independent, we can write 

±A(n) = Q(n) ± iU (n) 

= ^ ' ±Xim ±2 I^)m(fi) 

lm 

= YXE lm ± iBim) ± 2 Yim(n) , (42) 

lm 

where ± 2 h)m(n) are spin-2 spherical harmonics [14, 15]. 
af m = Ei m and af m = Bi m are the E-mode and B-modes 
spin spherical harmonic coefficients. The Stokes parame¬ 
ter Q and U are real and hence the complex conjugate of 
the spin-2 field reduces to + X(n)* =_ X(n). The spin 
spherical harmonic coefficients can be expressed as 

jj, +Xim -f-__ Xim 

-t^lm = ^ 

m --Xi m / ao\ 

Bim = -—- ■ (44j 

The complex conjugates of the above coefficients are 

+ x *m = (-1 rX^ m 
EL = (-1 ) m E lt _ m 
BL = (-1 ) m B lt _ m . (44) 

We briefly layout the formalism of polarization BipoSH 
coefficients in the next section. In the following sections 
we discuss only the results of the first order lensing cor¬ 
rections to these coefficients and their contribution to the 
polarization angular power spectrum. 


C(n,n') = Y i A hhh Y h( A ) ® 2 Y h (n')} L M 

LMhh 

+ 2Af^{- 2 Yi 1 (n) ® _ 2 Yi 2 (n-')}LM 
+ 3 A iit*{2Y h (n) ® - 2 Yi 2 (n')}LM 
+ 4 Af^{- 2 Y h (n) ® 2 Y l2 (n')} LM ) , (46) 

where are the BipoSH coefficients for the corre¬ 

sponding BipoSH basis. Comparing Eqs. 45 and 46 and 
parity properties of the coefficients as in Eq. 44, it can 
easily be shown that 1 Af'J* = 2 Af^ and = 4 Afjf. 

The + and the - signs on the spin-BipoSH coefficients 
denote the angular averages of like and unlike terms, viz. 
((+X+X), {-X-X)) and ((+A_A), (-X+X)) respec¬ 
tively. These coefficients can be expressed in terms of the 
spin-2 spherical harmonic coefficients in the same way as 
for the temperature case. 


+ A hl 2 — y 1 {(+A; lTOl +A; 2 T 7 j 2 ) + (—A/ imi -A; 2m2 )} 

m±m2 

C LM 

^ '-'lim 1 l 2 m 2 


— y ] (( Fi imi Ei 2rn2 Bi imi Bi 2m2 )Ci irni i 2m2 


rpLM tdLM 


(47) 


where EPf* and BPf* are the E-mode and B-mode Bi- 

lll2 _ tlt 2 

poSH coefficients respectively. Similarly, 


-A-lih ~ {(+^1”1 -Xl 2 m 2 )(—Xi imi + Xi 2 m 2 )} 


C L 


LM 

m\l2m2 


rpLM 


+ B; 


LM 

Ul2 


(48) 


Under the assumption of SI, the spin-2 BipoSH coeffi¬ 
cients reduce to combinations of the standard E-mode 
and B-mode angular power spectrum as 


+<q = (—l^H^-CT) 

- A hh = (-l^H LCY + C™). (49) 
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VIII. DISCUSSIONS 

The BipoSH basis forms a natural basis to study the 
CMB two point correlation function of the CMB sky. 
The angular power spectrum is a subset of the BipoSH 
spectra. While a SI Gaussian CMB sky is fully charac¬ 
terized by the angular power spectrum, a non-SI CMB 
sky is characterized by the whole set of BipoSH spectra. 
In a SI universe, weak lensing due to LSS surrounding 
the observer distorts the CMB sky and induces modifica¬ 
tion to the CMB angular power spectrum which can be 
evaluated given the angular power spectrum of the pro¬ 
jected lensing potential ip. In this article we study lensing 
in a non-SI universe, in which the statistics of the CMB 
sky and the projected lensing potential are specified com¬ 
pletely by the BipoSH spectra. 

Analogous to the lensing modifications to the CMB 
angular power spectrum in a SI universe, we evaluate the 
lensing modification to the BipoSH spectra due to weak 
lensing in an non-SI universe. We provide an equation for 
the BipoSH coefficients of the lensed CMB sky evaluated 
in terms of the BipoSH spectra of the unlensed CMB sky 
and those of the projected lensing potential. We study 
the result under various approximations. Firstly we note 
that in the case of SI CMB sky, lensing of the CMB by 
the particular realization of LSS surrounding the observer 
generates BipoSH spectra. We further show that the 
bipolar map and the bipolar power spectrum proposed as 
blind tests of SI violation relate to the projected lensing 
potential and it power spectrum respectively. We note 
that these estimator are closely related to the optimal 
lensing estimators [11]. 

Next we evaluate the corrections to the CMB angular 
power spectrum due to the coupling between the non-SI 
power in the projected lensing potential and the primor¬ 
dial CMB sky. It is found that this non-standard cor¬ 
rection vanishes if either the projected lensing potential 
or the CMB at the surface of last scattering is SI. This 
implies that any correction to the power spectrum is at 
second order in the parameter characterizing the isotropy 


violation. This analysis also leads to conclusion that any 
late time isotropy violating phenomena e.g direction de¬ 
pendent expansion, which only makes the lensing map 
non-SI cannot result in any corrections to the angular 
power spectrum of the observed sky. 

We also provide an expression describing these non¬ 
standard corrections to the angular power spectrum in a 
non-SI universe. A search for these non-standard correc¬ 
tions in the measured power spectrum can in principle be 
used to indicate deviations from SI. For instance, the ob¬ 
served temperature angular power spectrum from recent 
measurement has considerable deviations from standard 
ACDM angular power spectrum at certain multipole mo¬ 
ments like l = 2, 22 and 40 which could be potentially 
indication of deviations from SI. 

We also study the effects of the higher order lensing 
correction terms in a SI universe. While these correction 
terms does not contribute at the angular power spectrum 
level in a SI universe, they introduce some additional 
terms in the variance. We evaluate this correction to 
the cosmic variance in a SI universe due to weak lensing 
effects. 

Finally we discuss the BipoSH formalism for describing 
CMB polarization in a non-SI universe. The polarization 
BipoSH basis consists of combinations of spin-Spherical 
Harmonics and the corresponding BipoSH coefficients are 
further characterized as E-mode and B-mode BipoSH co¬ 
efficients. We can further calculate how these coefficients 
are modified due to weak lensing as in the temperature 
case and it can be shown that the modifications to the po¬ 
larization BipoSH coefficients closely resemble the mod¬ 
ifications to the temperature BipoSH coefficients. 
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Appendix A: Calculation of correction terms to temperature 


BipoSH coefficients due to lensing 


The two point correlation function for the lensed temperature field can be expressed in terms of the unlensed 
temperature field and the projected lensing potential through Equation 20. Each of the terms in the two point 
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correlation function of the lensed temperature field in Equation 20, can be expressed in the harmonic space as 
Y A hhi Y ‘A^i) 0 Yi 2 (h 2 )} L M= E A hh{ Y h( A i)®Y l2 (n 2 )} L M 


LM 

hh 


LMhl 2 


+ J2^’ knY kn( fi l) E A hh'{ Y h:a(ni)®Y h (ii2)}L<M' 

kn L'M'lsh 

+ ^^kvXkni fi 2 ) J2 A hh'{ Y h( A l) ®Y U :a (n 2 )}z/M' 


kn 


L'M'hh 


+ E ^k^iYki (ni) ® Y k 2 {n 2 )}KN E A hh'i Y h-A ni) 0 Yi 4 :b (n 2 )} L , M ' 


KN 

k\k2 


L'M' 

Z3Z4 


+ 9 ^ ^f , {^3:a6(ni)0^ 4 (n2)}L'M' 


KN 

k\k2 


L'M' 

/3Z4 


+ 7 ; E (* 2 ) 0 >fe 2 (ni)}icjv X! A hh' { Y i3--ab(n 2 ) 0 V; 4 (ni)}z,'M' (Al) 


KN 

k\k2 


L'M' 

Z3Z4 


Using the orthonormality relation of the BipoSH 

JJ dn 1 dn 2 {Yi 1 (n 1 ) 0 Yi 2 (n 2 )}LM{Ai 3 (ni) 0 Y u (n 2 )}* L , M , = Si^S^iJ LL '6 M M', 


(A2) 


the lensed temperature BipoSH coefficients NY can be obtained in terms of the BipoSH coefficients for the unlensed 
temperature field Aj' 3 u an d the projected lcnsing potential as 


~\LM aLM I „,LM { / \L'M’\ . „,LM / / \L'M'\ , ntm /,-v.K.N aL,m\ 

A hh — A hh + [ l a hh Vpkn, A lll2 )+ 2 <^hl 2 Nlm, A llh ) + Pzp 2 N k 3 k 2 > A ; 3 ; 4 ) 

LMr.r.KN aL'M'\ , LM /,t,KN \L'M'\ 

l7q/ 2 Aj 3 ; 4 )+ 17l'bfc 1 fc 2 i A /3 j 4 ) 


(A3) 


1. Calculation of the correction terms Ctr',; 

The correction terms Qq and Ct 2 are linear in the lensing deflection field and are expressed as 


= JJ dn ± dh 2 V“V’(n 1 )V^ 1 (T(n 1 )T(n 2 )){y Zl (n 1 )0F /2 (fi 2 )}2 M 
= [[ dnidn 2 Y^nYkni fi l) E E C W 4 m 4 y i3m 3 :a(ni)ii 4 m 4 (n 2 ) 

J J kn m 3 m 4 

X E Ci 2T n 2 ^ mi (fil)^(^) 

77117722 

= E^™ E E jd&i Y£Y hrna .. a YC imi 

kn L'M'IsIa 77117712 m 3 7714 

x J dn 2 Y hmi Y* m2 

= E^ E <*5 E ^'U/mil‘ 2 m‘N'l 3 , m 3 i 2 m 2 J ^nb| 3m3;a lj( mi 

kn L'M'lz m\m 2 mz 


(A4) 


= E^» E ^ rminm3 


/cn L'M'h 


mim2m3 


L M nLM ttt, 

^'I^m^l2m2^'lim\l2m2 1 l\ 


(AS) 
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Similarly 


n LM _ \ ' I \ aL’M' \ pL'M' riLM rm 2 ram 3 

2U hh — Z^^ kn hh L 'hm 3 hm 1 L 'hm 1 l 2 m 2 1 l 2 kl 3 


(A 6 ) 


kn L'M'l 3 


n/LM _ LM | r^LM 
a hl 2 - l a hh + 2 Ot hh 


(A7) 


2. Calculation of the correction terms f }j ; 

The second order correction term is quadratic in the lensing deflection field and is expressed in terms of the 
BipoSH coefficients of the lensing deflection field, 


PlZ = II [V^V" 2 (T(n 1 )T(fi 2 ))] {Y h (n,) ® Y h (h 2 )}* LM 

dnid n 2 V? X V | 2 £ *££ £ 

KNk 1 k 2 n 1 n 2 


£ A&f CliXumJh^^YumA^) 

L'M'I^Ia 77137714 

E 


E Ct^mY^n i)^ m 2 (n 2 ) 

m\m 2 


E aL'M' \ ^ nLM nL'M' jm\n\m 3 jm 2 n 2 mA /a o\ 

^hh 7 j ^ l\m\l 2 m 2 '-■'Z37713Z47714 1 l 1 k\l 3 1 l 2 k 2 lA ' 

KNkik2 n\n 2 L'M'l 3/4 77117712 7713 7714 


3. Calculation of the correction terms 7 /^ 

The two components of the other second order correction term 7 ^f 2 7 are also expressed in terms of the BipoSH 
coefficients of the lensing deflection field, 


..LM 

1 7z 1 i 2 


Similarly, 


d&i dn 2 [V^V^^CAr^^i))] [V" 1 V" 1 (T(n 1 )T(n 2 ))] {^(n,) ® Y l2 (n 2 )}* LM 
dhi dh 2 V a Al V b Al J2 Kk2T / Ck 7 1 k2n 2 Y kl n 1 (^l)Y k2 n 2 (Al) 

KNk lk2 ni" 2 


V^V" 1 A hlf E ^XumJl3m3^l)Y limi {^) 

L'M'IzIa 7713 7714 

E ^iCAT aL'M' V' C LM nL' M' 

^ k\k 2 7 v ^'kinik 2 n 2 / v ^hh 7 ^ ^ l\m\l 2 m 2 Z37713Z47714 




KNklk2 711712 


L'M'hU 


77117712 7713 7714 


x J dniyfc“ ni (ni)y^ na (ni)yj 3m 3 :ab (ni)Y ; * roi (rh) J dn 2 Y Umi {h 2 )Y* 2m2 (n 2 ) 


\ ' vj/tfAT V- nKN ST aL'M' ST' siLM sjL'M 1 jm x n 

n 7 , v fcife 2 / y ( -'/cinifc 2 ri2 / y ^hh 7 i WimiZ 2 m 2 '-'Z3m 3 i 2 m 2 t -'i 1 fciZ 1 


7117713712 
&2 


KNk\k2 711712 


L'M'h 


(A9) 


LM _ \ ^ .lyATAf \ ' nKN \ ^ aL'M' \ A zjL 

27/iZ 2 — o / ^ ^feife 2 / > ^Zei7ii/c 2 7i2 / v ^Z 3 Zi / > Wi 


LM rL'M' 


7711/27712^/37713/17711 u l 2 k\l 3 k 2 


jm,2r 


KNk\k 2 711712 


L'M'h 


(A10) 


LM _ LM , LM 

7iii 2 - l7hl 2 + 27i li2 ■ 


(All) 
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Appendix B: Evaluating the integrals 


The integral appearing in equation 24 can be evaluated by performing integration by parts and the using the identity 

V^ m (n)i-Kr + l)^;(n) q asin[16], 


lZ?r 2 = f dnYr m (n)Y^ mi (n)Y hm2 .. a ( n)) 

= 7}[h(h + l) + h(h + l) ~ + J dhY^ m {h)Y hmi (h)Yi 2m2 (h)) 

= i [h(h + 1) + h{l 2 + 1) - 1(1 + 1)] J^cl% 20 Cl™ mihm2 . 


(Bl) 


The second order correction term 7 ^]^, contains an integral as e< T ua ti° n 28 where 


jm 1 n 1 m 3 n 2 
^ l\k\l 3 k 2 


= J dn 1 ^* mi (ni)y^ ni (ni)y A ; 2 6 ri2 (ni)Yi 3 m 3 :a &(ni). 


(B2) 


This integral can be expressed in terms of spin spherical harmonics s Yim. We choose a convenient basis to study spin 
spherical harmonics 


m± = — (eg =F iej,). 


(B3) 


It is seen that m+ and ro_ are complex conjugate of each other i.e. ml j_ = m - F and act as the lowering and raising 
operators respectively [14], 


m + • V s Ylm — sj { +S)( 2 + } s-lYlm 

n v _ /(l — s)(l + s + 1 ) v 
10 _ ■ V sYlm — y 2 s+1 Yi m . 

In this basis, the derivative of a spherical harmonic is given as 

Vy Zm = \J 1 ^ 1 + ^ [m + tYirn - m- {Yi m } . 

Using these notations, the ,7-integral can be expressed as 


(B4) 


(B5) 


jm\n\m 3 n 2 

J hk 1 l 3 k 2 




Q*ki ^k 2 J 


dn Y * imi I m +a l Y kini - m _ a _1 y feini ] [m +b i Y k2U2 - m_ 6 _i Y k2 n 2 ] ^3777.3 :ab 
dn Yi imi [ i Yk lUl iYk 2 n 2 T —iYk ini — 1 Y]$ 2n2 m_ a m_^ 

iYk 2 n 2 m_|_ a iii_5 —\Y} Z ini llfc 2 n 2 m -a m +6] Y l 3 m 3 :ab , 


(B6) 


noting that a*, = \Jk(k + l)/2 . To evaluate this integral, we utilize identities for spin spherical harmonics from [10] 
and [17] 


iUfc ini (n) S2 Yk 2 n 2 (n) 


dn siYi irni S 2 Yi 2Tn2 S 3 Yi 3m3 


E B-fcifc 2 rtk—s sikn 

rj— j-r U k 1 -s 1 k 2 - S2 U k 1 n 1 k 2 n 2 
kns * k 


sYkn (H) 


n; 2 n; 3 


cf i_Si , c 

12—S213—S3 


Zl?7li 

^2 777.2^3 777-3 ' 


(B7) 












13 


The first term of the j7-integral simplifies to 

\J\ g fci g fc2 j dfk ^mi l^fcini l^k 2 n 2 iV^-+a^ )(^+Z>^ )Uz3m 3 

= a kl a k2 J2 -^r^r-ifa-iCn^n, / dfi ^ A(m+»V“) x (-a, 8 ) 

lb. 7 

= ~«fci«fc 2 a l 3 El yZ-j^n C fc 1 -lfc2-A~imfc2n2 


kn 

_ I \ ' n fc!fc 2 r>k —2 pkn ^-,/ 1 Q /i^mi 

Uk 1 U‘k 2 U'l 3 Ul 3 2 _^ rj- „ L 'fc 1 -lfc2-l L 'fcinife 2 ra2 rr~ j-r k- 2 l 3 2 ^knl 3 m 3 

kn V k V Ii 


Q j k\Q j k 2 Q j i 3 bi 3 n fcl n fe 2 n i3 v ' ^, k ~2 r>i\ 


E /ik —2 /-jZiO \ ' v /'jfcn /^Zi 

L 'fci-lfc 2 -l U /c- 2 Z 3 2 L, feinifc 2 n 2 u fer 


mi 

/cnZ 3 m 3 ’ 


Since the Clebsch Gordon coefficients j are valid only for m = m\ + m 2 , the sum over s gets fixed to 
in this case. Similarly the other terms in the J integral are 


J-i = 


g fcl g fc 2 g l3^3 Bfc 2 Bz 3 ^ ' ^k 2 y^liO \ ' nkn 
~ jj / J ^k 1 lk 2 l'~'k 2 l 3 -2 / ^k 1 nik 2 n 2 L 'knl 3 m 3 


rr _ ^ZiO nkn n h 

J3 ~ 47r / , u fci-lfc 2 l u fc 0 l 30 / , 4 irnbn 3 L b 


mi 

AmZ 3 m 3 


_ g fci g fc 2 a f 3 ^ rk0 i l0 rkn r i x 

^ ^ u fcilfc 2 — l u fcQZ 3 Q / v l/ fcinifc 2 n 2 u /cr 

1 k n 


m 1 

ZcnZ 3 m 3 


Using the identity from [12] 

nlm _ / 1 \Zi +Z 2 —Z nl—m 

^ZimiZ 2 m 2 V - 1 / ^Zi—miZ 2 —m 2 ’ 

the ^/-integral in equation B14 simplifies to 

&k\ ^fc 2 ^z 3 n fcl n fc 2 n z 


jm 1 n 1 m3n 2 _ _ 

J l\k\l3k2 


Hz, 


E hCi fe 2 -4oz30 {! + (-l) fcl+fe2 - fc } 


_i_ /’fcS ^ZiO r 1 I / 1 'ifei+fc 2 +l 3 -ll 'll \ ' fkn 

~r D h i 'k 1 lk 2 l u k 2 l 3 -2 l 1 + l _i l / / TtimfaB, L h 


imi 

■ / fcnZ 3 m 3 


If the unlensed CMB temperature field and the lensing deflection field is assumed to be SI, i.e. L' = M' = 
K = N = 0, the second order correction term 7 in equation 27b reduces to 


-,LM _ 1 \ ' ( 1 \ni \ ' pLM \( 1 \m 2 TTn 1 n 1 -m 2 -n 1 , / nmiri rm 2 »i-n 

1I1I2 9 / E ^Zci / ^ ^ZimiZ 2 m 2 [A ^ likil 2 ki 'V 4 / l 2 kil\ki 


kini m\m 2 

where we need to evaluate the quantity, 


E / 1 \ni T-mini -77 

l -W J ZifciZ 2 fci 


tt/ci«fci«Z 2 nfciiifejnzj. 


47T 


Hi! 


E 


the J integral 


T _ \ _J^kik2_„ k — s 

J 1 u fci u fc 2 / , /U”TT U k 1 -s 1 k 2 -s 2 '-'k 1 nik 2 n 2 

kns v ^ 


Qkn 


x J dn yj* mi s Yfcn [m +Q m +f) + m „m_ f( - m +a m_z, - m_ a m +6 ] ^Z 3 m 3 :ab 

life 


E 1 Zci/c 2 /7/cO /7fcn 


kn 


V^rllfc 


■'/ci0fc 2 0 v ^fcini/c 2 n 2 


1 . 


(B 8 ) 
s = 2 

(B9) 

(BIO) 

(Bll) 
0 and 

(B12) 

(B13) 


^fci ®fc: 


(B14) 
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The first term in the integral X is 

h = I dn Y* m Y kn (m +Q V°)(m +b V% m3 ) 

= -ai 3 j" dh Yi imi Yk n (m+ a V -iVz 3 m 3 ) 

(l 3 -l)(l 3 + 2) J dA Y * m y kn _ 2 Y i3m3 = 0 (B15) 

The above term vanishes for the integral of product of two spin-zero and one spin 2 spherical harmonic functions. 
The second term in the integral X vanish by similar arguments. The third term in the integral turns out to be same 
as the fourth term and is given as 



% = - j dnY* m Y kn (m +a V Q )(rn_ b V% m3 ) 
= ~ a ^ 3 j^ kn (m_|_ a V 1 l; 3 m 3 ) 

= 0; 3 J dn Yi imi Ykn Yi 3 m 3 = I 4 . 

Using the above equation and the identity from [12] 


dn Y* V Y, — i° rhm! 

U11 1 hm 1 I hm 2 I hm 3 — ^hOlsO^hmilsms’ 


\/47rII ;i 


the j7-integral in equation B14 simplifies to 

®;rr 2 = ^^ a 3 y . 


2 \ 1 nfc 1 IIfc 2 II; 3 k Q rtli 0 rilimi 

a tt ^'ki0k20'-'kinik2n2'-'k0l 3 0^'knl 3 m 3 ' 

kn ll 


(B16) 


(B17) 


(B18) 


Appendix C: Calculation of the correction terms in the cosmic variance 


The variance of the lensed CMB angular power spectrum to first order terms in Cf is given by 


c t 


(C 1 C 1 ) - {Cif 

(cf) + + <A*°%°> + «°7S 0 )] + ^<< 0 <) 

(Cl) 2 + ^ m oo )(ffi°) + «°}<7H 0 )] , 


(Cl) 


where we have used the relation A™ = (—l) ; n/Cz and the correction terms a,/? and 7 are given in Equations 35a, 

~ 2 

35b and 35c. Firstly, we focus on the terms contributing to (C\) and find 


«°)(ffi) = (-1 ) l n l c l (-i) l ihJ2 c ii c > 


kh 


F(hki) n k n hrl0 

yfa n, k0h y 


1 2 


= c^CiA 


kh 


F(hkl) 

\[Yk 


n fc n ;i c“ 


1 2 


fcOZiO 


(A*? 0 ) (iii 0 ) = (-i) l n l c l (-i) l+ 1 n l i(i + i)c l R = -i(i + i)RCfuf. 


(C2) 


Next we evaluate the terms which contribute to (Cf) in Equation Cl. The term (A ; * ; 00 a™) reduces to zero as aff 
does not contribute at power spectrum level. The terms which contribute are the cross-terms (A*f a (iff). (A*® 0 7 $°) 
and the covariance term (|a™| 2 ). Using the covariance relation for bipolar coefficients from Equation 37, we find 

( A u°Kh) = Cf {d n Ai 2 + (-1 ) l+l+ 0 6 ii 2 6 i h ) S L 0 S M0 + CiCiAAA-iy^ShhdLo (C3) 
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Using the above Equation C3, the first cross-term is 


C-lVi—mi /_1 \l—m 

{ATfll) = ^- 1 ) nC t( 2C l S ^U 2 +C l C h U l U ll (-iy +h 6 hl ^ J2 ^- jjnrimi j—m—n—mi 


knli 


lb 


jj M lkl\ ± lkl\ 


= 2 


2 F(lkl) l0 

ll^C 


1 2 


\/47 T 


kOlO 


n? 


EE^^+E^^ 


kU 


Fihkl) 

\[Fk 


1 2 


n fc n ;i 4° 0il0 


i 

n 2 


= (A0 0 )(/3°°}+2E^’Q S 


F(lkl) rr ^,jo 


-\/47r 


■n fc c, 


fcO/O 


i 2 

/ v v V-'knlimi J 

m nmi 

(C4) 


where we have used the relation for sum of products of two Clebsch-Gordon coefficients from [12] 


n 2 

\ nhm3 fhm3 _ _ h_c c 

Vimi/2m2Uimi/4m4 jj 2 h^4^ m 2 m 4 

h 


m\rri3 


(C5) 


and the relation Y^ m = H?- The second cross-term can be evaluated as follows 


/ /I^OO-.OOx 

\ A u la ) = 


(-i y 


E n 2 wk 

knm L 

i * 00 \ /_,00 


/^VVo /^2 . yr 2 si 2 \ jmn—m—n 
U k V ZO Z “T ^ l )Ji 


J Iklk 


= -l(l+l){2Cf + CfHf}R 


= (ADm-2Cfl(l+l)R. 


(C 6 ) 


The covariance of the a term is given by 


K%r> 


A Tt 2 \ ' \ ' / [4 / \ /\kn A*k'n'\F(l\kl)F(^l-^kl) sjkO nk '0 

Z^ z^ n 2 n n 

knh k'n'l'i 1 k k ' 

= 4 E E M(<n + (-l) fl+i+fe ^i,i^t'i>] f "[n]; CIS, oC&jt 

4^C*Cf 


fcnZi k'n'l 2 


/c/l 


F{l\kl) Ilqllfc ;o 
V^Br II; 


fcoqo 


= 5 iK 0 ){A 0 1 0 ) + 4E^Q 2 


F(/fcZ) 

•\/47r 


fc 

n 2 


F(lkl) tt 


\/47r 


■n fc c 


n fc c i0 


fcO/O 


(C7) 


Note that the cosmic variance for the unlensed CMB power spectrum is given by (C 2 ) — (C ;) 2 = 2C?/(2l + 1). Also 
the lensed CMB power spectrum is 


Ci = Ci - Rl(l + l)Ci + E G k C h F{ ' h \ k ' l) 

4?r 


4Ifcb[i pi o 

n* 


n 2 


■'fcoqo 


(C 8 ) 


Retaining terms which are first order in the lensing potential 


Ci « Q 2 -2MG + l)Cf + 2C' i EC'fcC'b 

kl\ 


F(h,k,l ) 2 

47 T 


n fc n, ^, ;o 
n. 


■'kOhO 


(C9) 
































16 


Using the equations C2, C4, C 6 and C7, the cosmic variance for the lensed CMB power spectrum for terms which are 
linear in Cf can be expressed as follows, 


Ci 


- ( C f) + 772 
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«)(/3°°)+2^C^ 2 
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- 2 Cfl(l + 1 )R + ^(A™)^ 0 ) + 2 Cf Y C, 
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F(hkl) IUII; l0 
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